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摘要 ” 托 卡 马克 等 离子 体 破裂 会 产生 逃逸 电流 , 如 不 进行 抑制 , 其 携带 的 巨大 能 量 将 对 设备 造成 严重 破坏 。 


本 文 使 用 DREAM 程序 中 的 流体 模型 ， 研 究 注 入 掘 毛 / 氛 混合 气体 对 破裂 逃逸 电流 的 影响 。 研 究 表明 ， 注 入 


气 氢 / 氛 混合 气体 可 以 抑制 最 终 形成 的 平台 逃逸 电流 。 在 讨论 的 破裂 前 等 离子 体 电流 p YEH 


内 ， 最 优 条 件 下 
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氮 / 氛 在 混合 气体 中 的 含量 应 在 0.50% ~ 0.70%， 气 的 注入 量 应 在 1029 m3 ~ 10?! m3. ZEXXA GB 
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效果 都 会 减弱 ， 甚 至 会 增 大 逃逸 电流 。 破 裂 前 等 离子 体 电 流 Ip 是 影响 逃 
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级 托 卡 马克 装置 上 , 注入 混合 气体 的 密度 需要 达到 10? mo, 这 是 目前 大 量 气体 注入 (Massive Gas Injection, 


MGI) 技术 所 不 能 达到 的 ， 通 过 散 裂 弹丸 注入 气 握 / 氛 混 合 物 将 是 更 加 可 行 的 方式 。 


关键 词 : ” 托 卡 马克 ， 等 离子 体 破 裂 ， 逃 逸 电流 ， 大 量 混合 气体 注入 ， 和 气 毛 /所 混合 气体 
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Abstract [Background] Tokamak plasma disruption generate runaway current that, if not suppressed, carry 
enormous amounts of energy that can cause severe damage to equipment. [Purpose] Investigate of the effect of 
injecting deuterium-argon/neon mixture gas on the runaway current during disruption. [Methods] In this paper, 
numerical simulations are carried out using a fluid model in the DREAM code, which is capable of self-consistently 
simulating the evolution of plasma parameters with time during the disruption process. [Results] It is shown that 
injection of a deuterium-argon/neon mixture gas suppresses the eventual formation of platform runaway current, but 
the deuterium-argon/neon mixture gas usually has an optimal content and ratio. In the range of pre-disruption plasma 
current Jp discussed in this paper, the content of neon/argon in the mixture gas should be in the range of 0.50% ~ 
0.70% and the injection amount of deuterium should be in the range of 10?? m? ~ 107! m?. Outside this range, 
deuterium-argon/neon mixture gas injection has a diminished effect on the suppression of runaway current and even 
increases them. [Conclusions] The pre-disruption plasma current /p is the key factor influencing the runaway current. 
The larger Jp is the larger the runaway current is formed and more mixture gas needs to be injected. On fusion 
reactor-scale tokamak devices with Jp up to 10 MA, the amount of injected mixture gas needs to reach 10? m? 
which is not achievable with the current massive gas injection (MGI) technique and injection of a deuterium- 
argon/neon mixture through a shattered pellet would be a more viable approach. 
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argon/neon mixture gas 


利用 托 卡 马克 装置 产生 聚变 能 所 面临 的 主要 问题 之 一 是 等 离子 体 破 裂 。 等 离子 体 破 裂 会 
产生 “ 逃 揭 电子 ?0， 如 果 不 对 其 有 效 的 抑制 和 缓解 ， 其 携带 的 巨大 能 量 将 会 损伤 装置 内 壁 部 
件 与 设备 。 在 大 型 托 卡 马克 装置 ， 如 欧洲 联合 环流 器 (JET: Joint European Torus) ~ HÆ JT- 


60U 装置 (JT-60U: JAERITokamak-60 Upgrade) 、 德 国 轴 对 称 偏 滤器 实验 装置 (ASDEX-U: 


Axially Symmetric Divertor Experiment Upgrade) 和 扒 级 聚变 装置 ， 等 离子 体 电流 较 高 。 国 际 


热 核 聚变 实验 堆 CITER: International Thermonuclear Experimental Reactor) 的 等 离子 体 电 流 


可 达到 10 兆 安培 CMA: Megaampere) ， 破 裂 产生 的 逃逸 电流 高 达 几 MAP?I, 3r T TE 
放 的 高 额 能 量 会 破坏 装置 的 第 一 壁 ， 造 成 严重 的 损伤 ， 甚 至 会 缩短 装置 使 用 寿命 。 尺 管 在 目 
前 的 小 型 托 卡 马克 装置 上 ， 等 离子 体 电流 小 ， 约 1MA， 因 等 离子 体 破 裂 产生 的 能 量 和 热 负 
荷 对 装置 内 壁 的 损伤 可 忽略 不 计 外 , 但 在 中 小 型 托 卡 马克 装置 上 研究 探索 破裂 逃逸 电子 物理 
及 对 其 有 效 抑制 和 缓解 的 方法 仍然 意义 重大 ， 特 别 是 在 1 MA 以 上 的 大 等 离子 体 电流 条 件 
下 。 基 于 现行 装置 参数 开发 高 效 的 破裂 逃逸 电流 抑制 方法 将 会 应 用 到 未 来 更 高 的 等 离子 体 电 
流 托 卡 马 死 聚变 装置 上 。 

在 托 卡 马克 等 离子 体 破裂 过 程 中 , 等 离子 体 的 热能 迅速 散失 ,导致 热 狸 灭 (TQ: Thermal 
Quench) 现象 的 出 现 ， 随 着 等 离子 体 的 冷却 ， 电 阻 急 剧 上 升 ， 从 而 引起 等 离子 体 电流 快速 下 
VEJÉ FAILS. CCQ: Current Quench) 。 在 电流 狂 灭 的 过 程 中 ， 逃 逸 电子 通过 初级 机 制 产 
增值 。 初 级 机 制 主要 包括 Dreicer 机 制 中 5 、 氛 衰变 产生 的 1 射线 、 活 化 壁 
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发 射 的 高 能 ”射线 的 康 普 顿 散射 和 热 尾 机 4 


逸 种 子 通过 次 级 雪 衣 机制 倍 增 形成 逃逸 电流 外。 


在 托 卡 马克 等 离子 体 中 ， 


的 电场 力 与 电子 的 运动 速度 无 关 ， 
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的 摩擦 阻 力 会 随 着 速度 的 增加 而 减 小 , 但 电子 受到 


因此 当 电 子 速度 增加 到 一 定 程度 时 , 电子 受到 的 阻力 将 会 


小 于 受到 的 电场 力 , 这 使 得 电子 被 持续 加 速成 为 逃逸 电子 。 当 等 离子 体 破裂 发 生 时 ,快速 衰 
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减 的 电流 会 引起 感应 电场 已 的 产生 , 这 个 电场 通常 远大 于 子 临界 电场 五 = 7 DIL 
Z£m,c 


中 ，e 


程 中 ， 


为 电子 电量 ，InA 为 库伦 对 数 ， 它 是 等 离子 体温 度 和 密度 的 缓 变 函 数 ， 对 于 磁 约 束 等 离 
子 体 一 般 取 16~17，ne 为 电子 密度 ，s0 为 真空 介 电 常数 ，me 为 电子 静止 质量 ，c 为 真空 中 的 
光速 。 热 电子 在 与 逃逸 电子 的 碰撞 过 程 中 获得 能 量 ， 使 得 更 多 电子 成 为 逃逸 电子 。 在 这 个 过 


逃逸 电子 数目 以 指数 型 快速 增长 ， 这 就 是 逃逸 电子 的 雪崩 倍 增 过 程 。 


为 了 缓解 等 离子 体 破裂 @8 产 生 的 危害 , 目前 常用 的 方法 有 


Ti: 大 量 气 体 注入 (MGI: 


Massive Gas Injection)、 杀 手 弹 丸 注入 (KPTI: Killer Pellet Injection ) 以 及 散 裂 弹丸 注入 (SPILD: 


Shattered Pellet Injection) 。 这 些 方法 通过 向 等 离子 体 中 注入 气体 或 小 颗粒 来 消耗 逃逸 电子 能 


Hi 


TATAG. MGI 是 目 
Experimental Advanced Superconducting Tokamak), "P 
U 等 装置 上 者 
气体 往往 很 难 有 效 缓解 破裂 后 的 逃逸 电流 ， 通 过 加 入 少量 的 高 原 


， 增 强 等 离子 体 的 稳定 性 和 抗 干扰 性 能 。SPI 被 应 月 


以 更 力 


HI 有 效 的 抑制 逃逸 电流 。 数 值 模 执 胡 
12 MA 的 逃逸 电流 ， 而 月 


FO; MGI 设施 ， 可 用 于 抑制 破裂 产生 的 逃逸 电流 。 


在 少 部 分 装置 ， 但 SPI 设计 和 优化 具 


有 挑战 性 ， 较 大 的 固体 弹丸 可 能 会 损伤 装置 第 一 壁 ， 需要 精心 设计 , 保证 其 大 部 分 杂质 能 在 
前 普遍 采取 的 方法 ， 在 全 超 导 托 卡 马 克 东 方 超 环 (2 (EAST: 


国 环流 器 二 号 A (HL-2A) 和 ASDEX- 


但 是 ， 通过 MGI 注入 单一 


子 序数 (高 Z) 杂质 气体 可 


究 发 现在 ITER 参数 条 件 下 等 离子 体 破裂 会 产生 约 


日 MGI 注 入 气氛 混合 气体 后 ， 最 终 形 成 约 3.7 MA 的 逃逸 电流 [93]。 


因此 通过 MGI 注入 包含 一 定量 杂质 的 混合 气体 更 适合 于 抑制 破裂 逃逸 电流 。 


握 和 氛 是 常 被 用 于 注入 到 托 卡 马克 等 离子 体 中 的 杂质 气体 , 因 其 较 高 的 原子 序数 ,可 以 
高 效 地 通过 辐射 降低 等 离子 体 的 储 能 。 但 过 量 的 握 或 氛 会 导致 等 离子 体温 度 快速 下 降 , 更 容 


易 狂 发 热 狂 灭 ， 产 4 
好 的 与 燃料 相 容 ， 又 可 通过 增 大 等 离子 体 密 度 降 低 破 裂 逃逸 日 
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更 多 的 逃逸 电流 。 和 所 是 未 来 气氛 聚变 堆 燃料 之 一 ， 注 入 大 量 的 气 可 以 较 


捉 托 卡 马 克 等 离子 体 破裂 逃逸 


电子 的 能 量 。 因 此 ， 通 过 MGI 
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合 气体 各 自 的 含量 将 是 影响 抑制 效果 的 关键 因素 ,本文 基 于 中 国 环流 器 二 号 M CHL-2M) 托 
卡 马 克 装 置 大 等 离子 体 电流 放电 条 件 ， 研 究 大 量 气 毛 / 氛 混合 气体 注入 对 抑制 破裂 逃逸 电流 
的 影响 ， 并 对 比分 析 气 氢 和 气氛 两 种 MGI 对 抑制 破裂 逃逸 电流 的 效果 。 

本 文 将 使 用 托 卡 马克 破裂 逃逸 电子 模拟 程序 DREAM04 开 展 数 值 模拟 研究 , 该 程序 能 
自治 地 模拟 托 卡 马克 等 离子 体 破裂 过 程 中 等 离子 体 参 数 随 时 间 的 演化 ,在 第 一 部 分 介绍 本 文 
的 研究 背景 和 动机 后 ， 第 二 部 分 给 出 了 数值 模拟 方法 和 逃逸 电子 相关 的 物理 模型 。 第 三 部 分 
给 出 了 注入 不 同 含量 气 氮 / 氛 混 合 气体 后 对 逃逸 电流 的 抑制 效果 ， 并 进行 了 对 比分 析 。 最 


总 结 全 文 并 对 未 来 进行 展望 。 
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1 数值 模拟 方法 
等 离子 体 破裂 后 热能 会 在 几 毫 秒 之 内 释放 ,这 个 阶段 被 称 为 热 狐 灭 (CTQ)， 在 此 期 间 电 
了 温度 T 呈 指数 衰减 ， 温 度 从 初始 的 To(r) 向 最 终 温度 Ty(r) 衰减 ， 其 特征 衰减 时 间 常 数 为 


0， 温 度 的 指数 衰减 (2 是 由 方程 (1) 来 表示 的 


TG, t) -T,G) +T Cr) - T, Ge "^ a) 


ZRBSATACHS SSAA, A Seam, Mime CCQ) 。 
这 是 第 二 阶段 的 开始 , 电流 衰减 得 越 快 感应 电场 也 就 越 强 , 从 而 加 速 电子 使 之 成 为 逃逸 电流 。 
在 等 离子 体 中 电子 发 生 逃 逸 离 不 开 两 个 重要 的 电场 : 临界 电场 E Dreicer 临界 电场 


en, In A 


TO 是 速度 等 于 光速 的 电子 发 生 逃 逸 时 所 需 的 电场 。 当 温度 为 下 
4ze,m,c 


Ep. 临界 电场 ,= 


的 热电 子 受到 高 于 Dreicer 临界 电场 E, =en, In /4xe2T. 产 生 的 电场 力 时 ， 由 于 其 所 受 


到 的 摩擦 阻力 小 于 电场 力 ， 这 部 分 电子 就 会 被 持续 加 速 ， 从 而 发 生 逃 逸 05。 

等 离子 体 破 裂 后 有 许多 不 同 的 机 制 可 以 产生 逃逸 电子 ， 其 中 起 主要 作用 的 是 初级 
Dreicer 机 制 、 热 尾 hot-tail 机 制 以 及 次 级 雪 骨 机制。 在 这 里 我 们 忽略 了 热 尾 hot-tail 机 制 09， 
因为 在 本 文 研 究 中 ， 等 离子 体 在 破裂 后 有 足够 的 时 间 来 丧失 能 量 ， 冷 却 时 间 足 够 长 ， 整 体温 
度 下 降 , 高 能 尾部 的 粒子 会 逐渐 失去 能 量 , 因而 对 逃逸 电流 的 贡献 减弱 。 当 冷却 时 间 足 够 长 ， 
Dricer 机 制 会 占 主 导 地 位 ， 而 热 尾 机 制 在 这 种 情况 下 相对 较 弱 不 起 主导 作用 ， 所 以 忽略 了 热 
Fé hot-tail 机 制 。 因 此 ， 不 考虑 hot-tail 机 制 ， 总 的 逃逸 电子 的 密度 za 随时 间 演 化 的 方程 0 
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初级 机 制 中 ，Dreicer 机 制 产生 的 逃逸 种 子 07 
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这 里 的 z = 4eme | n,e* In A 为 相对 论 电子 碰撞 时 间 n8， 古 为 平行 电场 ， 逃 逸 种 子 通 


次 级 雪崩 机 制 倍增 09201 
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Z ， 其 中 nn BLZ, 分 别 为 离子 种 类 六 的 电荷 态 i 的 密度 和 电荷 数 。 这 里 


1 
0 = (1 + 1.46 €+ 1.72 €)7t, e-r/R. "ME Za l1, p= ARIES > Ec 时 ， 可 以 将 方程 
(4) 简化 为 : 


dn ELE —] 


( RE jo ~ (— 2 I 
m Tag 


dt 3rln A m 

Zar l 是 一 种 理想 情况 ， 通 常 Zer 会 大 于 1， 这 是 因为 杂质 (原子 序数 大 于 1) 的 进入 
会 使 得 Zor KT 1， 除 非 都 是 氧 同位 素 的 等 离子 体 。 实 验 上 ， 即 便 都 是 氨 同 位 素 等 离子 体 ， 

由 于 等 离子 体 与 壁 材料 相互 作用 产生 的 杂质 也 会 进入 等 离子 体 中 ,使 得 有 效 电 荷 大 于 1。 在 
本 文中 因为 注入 氛 氮 / 氛 混 合 气体 会 使 得 Zor 大 于 1, 所 以 , 在 本 文 计算 中 都 是 采用 方程 (4) 

进行 计算 。 

由 于 电子 、 气 离子 和 杂质 离子 的 碰撞 时 间 不 同 , 所 以 分 别 给 


出 这 三 类 粒子 的 能 量 平衡 方 


LE 


[20] 
30 (n,T) _ 3n, 7) p 
2 o CES 二 r2 + Pon — Par — Pion — Prine t 2: pt (6) 
3 Anaa) _ 3a ME prs poe : 


2 Ot 2r Or 


30 (njT,) 3n, ð OT, > oe 
ee eee ee ea, Pl ie i^ j 8 
2 Ot 2r Or (xr or )* p. dP: (8) 


其 中 Ne 为 电子 密度 ， X 


1 n? s! ZAP RAR, Py, = ojE? 是 欧姆 加 热 功率 密度 ， 


万 致 辐射 损失 功率 密度 Di ， 


rau 


P,, =1.69x 10% 6, p? D'JT, EV E, 


P, = >, BiTinin, 是 电离 能 量 损失 功率 密度 。 P? = 3n,(T, -T,)/ 2r; 为 能 量 平 衡 方程 中 


34247" e mm, T T; a 
不 同 种 类 粒子 之 间 的 碰撞 能 量 交 换 ,，T, = 77 A +—) yee Beary Te, H 
neZiZnA m m, 


中 下 标 i ，j 指 的 是 电子 以 及 所 和 杂质 离子 。 Py。 = 2 nz n, (n, T.) 是 线 辐射 功率 密度 ， 


它 是 所 有 电荷 态 密度 的 总 和 ， 其 中 辐射 率 工 。 从 ADAS 数据 库 中 获得 31。 同时， 通过 每 种 离 


子 种 类 的 每 种 电荷 态 的 密度 随时 间 演 化 的 方程 来 计算 杂质 的 电离 率 [17 


k 
dn, 


TUR n, nz. -(I, +R, )n;, 4 Ryan") (9) 


EB 为 第 个 电荷 态 的 电子 碰撞 电离 率 ，Ri 是 辐射 复合 率 。 


be 
1m 


随 着 等 离子 体 的 冷却 和 其 电导 率 的 下 降 (oj - 777), AV SB PS BPRS 


Ht 


子 体 电流 。 其 中 电流 密度 分 布 随时 间 的 演化 和 电场 的 演化 由 方程 (10) PARR 


CRAT 
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其 中 , k 为 拉 长 比 ， 等 离子 体 电流 密度 j = oS EJ, e o, E, ecn,, 为 欧姆 电流 密度 和 逃 


逸 电流 密度 之 和 。 因 是 逃逸 电流 密度 ，cw =3s0(2z7) /(0.51m;"e In A) Æ Spitzer Hi 


FR, Ly Bl c 分 别 为 磁 导 率 和 真空 中 的 光速 。 


基于 上 述 物理 模型 ,本 文 利 用 数值 模拟 的 方法 , 在 模拟 程序 中 设 定 总 模拟 时 间 和 最 大 和 迭 
代 次 数 , 并 当 满 足 任 一 条 件 时 , 方程 组 的 演化 将 终止 , 研究 MGI 技术 注入 气 和 高 Z 杂 质 〈 毛 
或 氛 ) 混 合 气 体 0222259 对 抑制 破裂 逃逸 电流 的 影响 。 在 程序 中 我 们 设 定 总 模拟 时 间 为 50 ms, 

在 这 个 时 候 结果 满足 模拟 参量 不 随时 间 演 化 。 当 注入 气体 的 密度 达到 102 m? 时 ， 碰 撞 耗 散 
效应 会 极 大 地 耗 散 逃逸 电子 的 动能 , 这 种 情况 下 不 会 形成 平台 逃逸 电流 ,而 是 形成 峰值 逃逸 


NS 


一 


外流 ， 即 逃逸 电流 达到 最 大 值 后 开始 快速 衰减 至 零 。 但 是 ， 现 在 的 MGI 技术 还 不 能 使 得 注 


入 气体 的 密度 达到 10? m? 量 级 。 因 此 ， 本 文 假设 注入 的 气 、 所 和 和 氛 的 密度 


范围 分 别 为 : np 


= 102 m3~ 10?! m3, n4; = 1018 m? ~ 2x10! m? 和 nye = 10!8 m3~ 2x10 m3， 并 且 注 入 的 混 


合 气体 是 均匀 的 。 


2 结果 与 讨论 


模拟 采用 HL-2M 托 卡 马克 装置 3 的 参数 , 装置 小 半径 qa=65.5 cm, 


大 半径 Ro= 178.1 


cm， 环 向 磁场 Bo= 1.8T， 破 裂 前 等 离子 体 电流 万 = 1.2MA。 初 始 电 子 温度 和 密度 的 径 问 


分 布 如 图 1 所 示 ， 中 心 电 子 的 温度 和 密度 分 别 为 Teo= 6.8 keV 和 neo = 4.6x10? m3。 破 裂 


后 等 离子 体温 度 快 速 下 降 ，TQ 期 间 的 电子 温度 衰减 根据 方程 (1) 演化， 参考 文献 [13]， 


W to=1ms, T= 10eV. 
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图 1 初始 的 电子 温度 (a) 和 密度 (b) 的 径 向 分 布 


0.8 1 


Fig.1 Initial radial distribution of electron temperature (a) and density (b) 


破裂 发 生 后 产生 逃逸 电流 , 并 形成 由 逃逸 电流 构成 的 逃逸 平台 , 图 2 


展示 了 在 破裂 后 没 


有 杂质 注入 时 1.2 MA 的 五 最 终 会 形成 约 0.8 MA 的 平台 逃逸 电流 ， 约 在 8ms 时 逃逸 电流 he 


与 等 离子 体 电流 及 相等 ， 即 逃逸 电流 完全 蔡 代 了 等 离子 体 电流 ,在 8ms 后 I. 5 I5 BEEN [8] 


变化 , 此 时 的 电流 定义 为 平台 逃逸 


外 流 应 。 和 破裂 所 带 来 的 危害 程度 在 一 定 程度 上 取决 于 平台 


逃逸 电流 与 五 的 比值 天/ 邦 。 在 图 2 中 ， 压 208 MA, IR) =0.75. XA 


克 装 置 上 是 不 能 被 接受 的 ， 如 ITER 装置 上 万 高 达 十 几 MA, 将 严重 影响 


比值 在 大 型 托 卡 马 
装置 内 壁 部 件 和 设 


备 的 安全 。 对 此 ,本 文 接 下 来 将 研究 在 破裂 时 注入 气 所 和 和 气氛 混合 气体 对 平台 逃逸 电流 的 影 


响 ， 给 出 有 效 抑制 平台 逃逸 电流 下 混合 气体 含量 和 平台 逃逸 电流 大 小 。 
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2 无 混合 气体 注入 时 等 离子 体 电流 和 逃逸 电流 随时 间 的 演化 


Fig.2 Evolution of plasma and runaway current with time without mixture gas injection 


2.1. 注入 气 氨 /所 混合 气体 的 结果 


为 了 理解 混合 气体 注入 对 抑制 逃逸 电流 的 效果 , 本 节 研 究 了 注入 气 握 或 气氛 混合 气体 对 


抑制 逃逸 电流 的 影响 。 在 破裂 后 注入 气 握 混合 气体 , 气 握 气体 的 含量 分 别 为 mp = 9.2x102m 


3 和 NAr = 1.1x10!9 m?, 图 3 (a) 展示 了 等 离子 体 电流 Ip. 欧姆 电流 Tonm 和 逃逸 电流 Tre 随时 
间 的 演化 。 可 以 看 到 注入 气 氮 混合 气体 后 形成 的 平台 逃 移 电流 显著 降低 ， 于 民 的 值 由 0.75 
降低 为 0.41， 这 说 明 混合 气体 的 注入 可 以 抑制 逃逸 电流 。 但 注入 此 混合 气体 含量 后 EL EG 


值 仍 相对 较 高 ， 因 此 需要 进一步 优化 注入 的 混合 气体 含量 。 图 3 CO 展示 了 破裂 后 电场 情 


况 ， 可 以 得 出 破裂 后 产生 的 电场 远大 于 E = 0.3 Vm， 此 时 中 心 电 子 密度 neo = 4.6x1019 m3. 


(85 č 
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3 YEATUMIRA AL (mp = 9.2x10% m?, na - 1.1x10 m?) 的 结果 。(a) 等 离子 体 电流 〈 实 线 ) ， 欧 


姆 电流 《点 线 ) 和 逃逸 电流 虚线 ) 的 演化 。(b) 不 同 归 一 化 径 向 半径 上 欧姆 电场 的 演化 。 


Fig.3 Results of injection of a deuterium-argon mixture gas ( np = 9.2 x 10% m?, nar = 1.1 x 10? m? ). (a) 
Evolution of plasma current (solid lines), ohmic current (dotted lines) and runaway current (dashed lines). (b) 


Evolution of the ohmic electric field on different normalised radial radii. 


El 4260 T^r 36r B TE A iad IR 4 Ca) 为 注入 气 的 量 在 np = 
9.2x10 m3 tf, 改变 氢气 体 的 含量 , 最 终 形成 的 平台 逃逸 
m3 时 会 达到 最 小 的 平台 逃逸 电流 , 下 =0.19MA。., 但 随 着 氢气 体 的 不 断 增 加 , 抑制 效果 减弱 ， 


平台 逃逸 电流 会 逐渐 增 大 。 图 4 b) 是 保持 na4x1075 m3, HUE UAE E Æ np=7x10” 


In 


电流 会 先 逐 渐 减 小 , TE na; = 4x1018 


m3 时 ， 平 台 逃 逸 电流 最 小 ,下 = 0.14MA， 随 着 气 气 体 不 断 增 加 ， 平 台 和 逃逸 电流 也 会 逐渐 增 


大 。 因 此 ， 可 以 得 到 最 优 的 气 毛 混合 气体 含量 : mp = 7x10? m? 和 mer= 4x10 m3。 在 此 优 


化 注入 条 件 下 ， 于 =0.14 MA 显著 低 于 没有 任何 抑制 措施 下 的 0.8 MA， 在 注入 的 混合 气体 


中 Ar 的 含量 nad/(npt+nar)x100% = 0.57%。 在 注入 的 气 握 混合 气体 中 ， 氛 的 量 过 大 会 导致 Ar 
的 含量 降低 ， 不 能 有 效 地 利用 Ar 来 辐射 损失 逃逸 电子 的 能 量 ， 但 是 过 量 的 Ar 会 进一步 降 
低 等 离子 体温 度 ， 促 使 更 高 的 感应 电场 产生 ， 从 而 产生 更 高 的 逃逸 电流 。 


np=9.2x1020m-3 | 


= 18-3 
n,,=4x10°"m 
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4 平台 逃逸 电流 随 注 入 气 握 含量 变化 的 结果 。(a) np = 9.2x1020 m3 不 变 ， 逐 渐 增 大 人 毛 的 含量 。(b) nar = 
4x108 m3 不 变 ， 逐 渐 增 大 气 的 含量 。 


Fig.4 Results of platform runaway current variation with injected deuterium-argon content. (a) np = 9.2 x 107° 


m? constant, gradually increasing argon content. (b) nar = 4 x 10!8 m? constant, gradually increasing the deuterium 


content. 
对 于 气氛 混合 气体 注入 ， 虽 然 其 对 逃逸 电子 的 作用 与 气 氨 混合 气体 一 样 ， 但 由 于 所 的 
原子 序数 比 氨 小 , 所 以 其 对 逃逸 电子 的 辐射 功率 损耗 相 比 更 小 。 图 5 中 展示 了 平台 和 逃逸 电流 
随 注入 气氛 含量 变化 的 结果 。 注 入 气氛 混合 气体 的 结果 与 图 4 中 注入 气 氮 的 结果 是 类 似 的 。 


在 最 优 条 件 下 ，np=7x102 m3, nye=4.5x10!8 m3， 平 台 逃 逸 电流 达到 最 小 值 ， 站 约 为 0.14 


N 


MA。 在 注入 的 混合 气体 中 Ne 的 含量 nne/(np+nne)x100% = 0.64%. 3B EG 5 S] 4 PRETT 


注入 的 量 , 气氛 注入 下 气 的 量 没有 发 生变 化 , 但 需要 注入 更 多 的 氛 , 或 者 混合 气体 中 氛 的 含 
量 要 稍微 增 大 。 
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5 平台 逃逸 电流 随 注入 气氛 含量 变化 的 结果 。 (a) np = 9.2x10 m3 不 变 ， 逐 渐 增 大 氛 的 含量 。(b) nxe 


RIT 


= 4.5x1018 m? 不 变 ， 逐 渐 增 大 和 气 的 含量 。 


Fig.5 Results of the variation of the platform runaway current with the injected deuterium-neon content. (a) np = 
9.2 x 10? m? constant, gradually increasing the amount of neon. (b) nye = 4.5 x 10!8 m? constant, gradually 


increasing the deuterium content. 


通过 对 上 述 结果 进行 分 析 , 可 以 得 到 注入 气 握 或 气氛 混合 气体 都 能 抑制 平台 逃逸 电流 ， 


但 同时 都 会 存在 一 最 优 混合 气体 含量 ， 当 超过 这 个 最 优 含量 之 后 ， 抑 制 效果 会 逐渐 减弱 。 由 


于 本 文中 注入 混合 气体 抑制 逃逸 电流 是 在 中 型 托 卡 马 殉 装 置 下 得 到 的 , 对 于 等 离子 体 电流 更 


大 、 和 破裂 前 等 离子 体温 度 更 高 的 堆 级 聚变 托 卡 马 死 装置 而 言 ， 可 以 忽略 装置 尺寸 的 影响 ,， 进 


步 研究 有 和 电子 温度 等 效应 的 影响 。 


2.2 ”破裂 前 等 离子 体 电流 和 电子 温度 效应 的 影响 


不 考虑 装置 尺寸 的 影响 ， 本 节 研 究 有 和 电子 温度 Tao 对 破裂 后 平台 逃逸 电流 的 影响 ， 并 


计算 给 出 不 同 克 和 Ta 情况 下 抑制 破裂 平台 逃逸 电流 效果 最 佳 的 气 氮 /所 混 合 气体 含量 。 
为 了 方便 对 结果 进行 对 比 , 表 1 给 出 了 增 大 五 或 Tio 时 在 没有 混合 气体 注入 的 情况 下 形 
成 的 平台 逃逸 电流 。 五 越 大 ， 最 后 形成 的 平台 逃逸 电流 六 也 越 高 ， 丰 你 随 着 五 的 增 大 而 增 


Ko MERK Teo 情况 下 ， 碟 稍微 减 小 ， 这 可 能 是 因为 过 高 的 Teo 使 得 等 离子 体 冷 却 时 间 尺 


度 延 长 所 引起 的 。 对 比 及 和 To 对 平台 逃逸 电流 的 影响 ， 万 是 影响 破裂 逃逸 电流 的 最 主要 的 


因素 。 这 对 于 五 高 达 ~ 10MA 的 聚变 托 卡 马克 反应 堆 而 言 ， 如 ITER， 必 须 采 取 有 效 的 措施 
来 抑制 破裂 后 形成 的 逃逸 电流 ， 减 小 其 危害 。 


注入 最 优 含量 的 气 氮 / 氛 混 合 气体 , 最 后 形成 的 平台 逃逸 电流 态 及 相应 V1 对 (1b, Teo) 


的 依赖 结果 列 于 表 2 和 表 3 中 。 从 表 2 中 可 以 发 现 , 在 注入 最 优 气 氢 混 合 气 体 后 , 平台 逃逸 


m 


EIIE f RKM. 1.2 MA 的 五 最 终 会 形成 约 0.14 MA 的 平台 逃逸 电流 ， 相 比 于 无 任何 


缓解 措施 下 的 0.8 MA， 逃 逸 电 流 在 气 握 混合 气体 注入 下 得 到 了 大 幅度 抑制 ，A&/4 的 值 也 比 


在 没有 注入 气 握 混合 气体 下 出 现 大 幅 下 降 。 但 是 , 当 石 增 大 一 倍 时 , 在 最 优 气 握 混合 气体 注 


AM elh WERTET KE. 而 Two 仍然 对 注入 气 握 混合 气体 后 的 平台 逃逸 电流 影响 


较 小 。 这 进一步 说 明 对 破裂 逃逸 电流 及 对 其 捉 秆 


— 


有 着 关键 的 影响 。 
Rl 无 混合 气体 注入 时 不 同 郊 和 Te 下 的 结果 


Tablel Results at different Jp and Teo without mixture gas injection 


Tp (MA) Teo (keV) IF. (MA) IP IT 
12 6.8 0.8 0.75 
1.8 6.8 1.4 0.76 
24 6.8 1.85 0.77 
1.2 10 0.6 0.5 
1.8 10 1.1 0.6 


EER MAMECK EEE p 变化 的 趋势 ， 对 于 有 效 抑制 


ITER 装置 等 离子 体 破 裂 后 的 逃逸 电流 ， 需 要 注入 氧气 体 的 含量 应 达到 102”m3， 毛 气体 含量 


lin 


的 范围 在 10? m? ~ 107° m? E BUE AA Bt ATE A CARAS ETE A 8 SEXA Bl] 10? m? 的 气体 。 


因此 , 对 于 未 来 的 高 等 离子 体 电 流 聚 变 托 卡 马 殉 反 应 堆 而 言 ， 需 要 更 加 优良 的 注入 方式 来 抑 


制 破裂 逃逸 电流 ， 比 如 散 裂 弹丸 注入 CSPI) 23, 


表 2 ”注入 最 优 含量 气 氨 混合 气体 的 结果 


Table2 Results of injecting the optimum content of deuterium-argon mixture gas 


nAy/(np+nar)x 100 
Ip (MA) Teo (keV) np (m?) nar (m?) 4 IÈ, (MA) Il 
(^j 


1.2 6.8 7x10? 4x10!8 0.57% 0.14 0.11 


1.8 6.8 1.4x10?! 8x10!8 0.57% 0.38 0.2 
2.4 6.8 1.9x10?! 9.5x10!* 0.5096 0.59 0.25 
1.2 10 4.5x1070 3.5x10!8 0.70% 0.12 0.1 
1.8 10 1.1x10?! 6.8x10!8 0.62% 0.31 0.14 
2.4 10 1.5x10?! 8.8x10!5 0.60% 0.5 0.2 


改变 混合 气体 中 杂质 的 成 分 ， 表 3 中 给 出 了 注入 气氛 混合 气体 时 的 结果 ， 结 果 表 明 对 
破裂 后 平台 逃逸 电流 的 抑制 效果 与 注入 所 毛 混合 气体 的 结果 相似 , 在 相同 气 注入 量 下 需要 注 
入 更 多 的 氛 才 能 使 得 最 后 的 ie 保持 基本 不 变 。 这 主要 是 因为 氛 的 原子 序数 比 毛 更 低 , 对 逃逸 
尼子 的 辐射 损失 能 力 更 弱 。 

X3 注入 最 优 含量 气氛 混合 气体 的 结果 


m 


Table3 Results of injecting the optimum content of deuterium-neon mixture gas 


NNe/(nD+NNe)X 10 


Ip (MA) Teo (keV) np (m?) nye (m?) IP, (MA) IE, /Ip 
‘0 
12 6.8 7x1020 4.5x10!8 0.64% 0.14 0.11 
1.8 6.8 1.4x1021 8.3x10!8 0.59% 0.37 0.2 
24 6.8 1.9x10?! 9.7x10!8 0.51% 0.57 0.24 


通过 对 上 述 结果 进行 分 析 和 研究 ， 可 以 得 出 在 中 小 型 托 卡 马克 装置 上 ， 其 五 普遍 在 1-3 
MA 以 下 ， 利 用 气 毛 / 氛 混 合 气体 注入 可 以 有 效 抑制 破裂 后 的 逃逸 电流 ， 且 注入 气体 的 含量 
不 会 超过 大 量 气 体 注入 技术 对 应 的 上 限 气 体 含量 。 在 优化 注入 条 件 下 , 从 表 2 和 表 3 中 可 以 
得 出 : 氮 或 氛 在 气 氮 / 氛 混 合 气 体 中 的 含量 应 在 0.50% ~ 0.70%， 气 的 注入 量 应 在 1079 m? ~ 
10?! m3， 这 些 可 以 使 得 逃逸 电流 最 为 有 效 地 被 抑制 。 在 最 优 注入 条 件 下 ， 被 抑制 的 逃逸 电 
流 仅 为 五 的 10%， 最 后 形成 约 100 -200 kA 的 逃逸 电流 。 因 此 ， 大 量 气体 注入 技术 ， 特 别 是 
气 握 / 氛 混合 气体 注入 对 于 抑制 现行 中 小 型 托 卡 马克 装置 上 的 破裂 逃逸 电流 是 有 效 的 。 但 是 ， 
IFMS p EE 10 MA 量 级 的 聚变 托 卡 马克 装置 上 的 破裂 逃逸 电流 ， 根 据 本 文 研 究 结 
进行 外 推 ， 需 要 注入 混合 气体 的 含量 应 达到 102? m3 的 量 级 ， 这 对 于 目前 大 量 气 体 注入 技术 
是 无 法 达到 的 。 


mi 
中 


由 于 大 量 气体 注入 技术 不 能 注入 含量 达 10? m? 量 级 的 气体 ， 对 逃逸 电子 的 碰撞 耗 散 效 


应 很 小 ,方程 (4) 中 EWE.c-1~ EWEce， 最 后 会 形成 平台 逃逸 电流 。 但 是 ， 如 果 能 够 注入 含量 


E3 


达到 102? m? 量 级 的 破裂 缓解 物质 ， 比 如 通过 SPI 技术 , 临界 电场 E 的 量 级 将 会 与 破裂 后 的 


欧姆 感应 电场 相当 ，EWE. ~ 1， 此 时 碰撞 耗 散 效应 将 会 对 破裂 逃逸 电流 产生 重大 影响 。 如 


图 6 所 示 ， 此 时 逃逸 电子 未 能 形成 平台 逃逸 电流 ,逃逸 电流 大 小 随时 间 先 增 大 后 迅速 减 小 至 


零 ， 在 约 19 ms 时 形成 峰值 逃逸 电流 ， 约 0.35 MA; 然后 在 约 70 ms 时 逃逸 电子 的 能 量 被 耗 


散 ， 逃 逸 电流 消失 。 如 果 考 虑 以 SPI 注 入 10? m? 量 级 的 氛 氮 / 氛 混 合 物 ， 则 形成 的 峰值 逃逸 


I 


m 


外 流 会 更 低 ， 而 且 会 更 快 地 将 逃逸 电子 的 能 量 耗 散 完 。 


(J 


对 于 托 卡 马克 破裂 而 言 ， 由 于 其 危害 巨大 ， 基 于 现行 装置 必须 发 展 有 效 的 预测 、 避 免 和 
抑制 技术 。 最 为 优先 的 应 为 预测 和 避免 技术 ,但 这 些 并 不 在 本 文 的 研究 范围 内 。 但 是 ， 即 使 


预测 和 避免 技术 来 防止 托 卡 马克 破裂 带 来 的 危害 , 发 展 有 效 的 破裂 抑制 缓解 技术 也 是 必要 


的 ， 因 为 预测 和 避免 这 些 技术 并 不 能 100% 地 实现 其 预定 的 目标 ， 破 裂 抑 制 技术 是 最 后 的 防 
御 措施 。 
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6 高 密度 (me = 107? m? ) 下 碰撞 耗 散 效应 对 逃逸 电流 的 影响 。 


Fig.6 Impact of collisional dissipation effects on runaway current at high density ( ne = 10? m? ). 


3 结语 


本 文 利 用 托 卡 马克 破裂 逃逸 电子 模拟 程序 DREAM 中 的 流体 模型 研究 大 量 气 氢 或 气氛 


混合 气体 注入 对 抑制 破裂 逃逸 电流 的 影响 。 通 过 对 比 有 /无 混合 气体 注入 条 件 下 平台 逃逸 电 


流 与 LAIR 的 大 小 ， 本 文 获得 了 和 气 


Hi 


Ha 


通过 


yey 
PI 


& BA TU A 


TE ROLE PRL p ARMA A 


MGI EREA ia TE TT EU IG 
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\ 体 中 的 最 优 含量 及 配 比 ， 并 在 最 优 
完 发 现在 最 优 注 入 条 件 下 被 抑制 的 逃逸 电流 仅 为 万 的 10%, 
制 当前 中 小 型 托 卡 马克 装置 破裂 后 的 逃逸 电 


等 离子 体 破裂 后 , 无 混合 气体 注入 会 产生 较 高 的 平台 逃逸 电流 , IPIS RK, 这 在 高 等 离 


子 体 电流 


KI ZA o 


量 及 配 比 之 外 的 范围 , 增 大 气 的 注入 量 会 导致 Ar/Ne TEI TE FA 
利用 ArNe 来 辐射 损失 逃逸 电子 的 能 量 ; 而 过 量 的 Ar/Ne 会 进一步 降低 等 离子 体温 


可 以 增 大 逃逸 电子 的 辐 

在 最 优 含 量 及 配 比 的 氛 氮 /所 混合 气体 注 
的 影响 最 大 ， 注 入 混合 气体 后 的 成 随 p 的 增 大 而 线性 增 大 。 注 
随 五 的 增 大 而 增 大 ， 当 五 = 2.4 MA 时 ， 注 入 混合 气体 的 总 量 达到 10?! m? 
围 内 ， 优 化 后 的 毛 或 氛 在 4 
以 此 外 推 


的 Tp 3 


at 


应 在 109 
合 气体 的 量 需 要 达到 10” m?, 

不 考虑 注入 方式 ， 本 文 看 
F 下 ,破裂 后 不 形成 平台 逃逸 


f 


在 此 条 伯 


的 大 型 托 卡 马克 装置 上 是 不 能 被 接受 的 。 不 


FRM 


EA AU lS nr UM n] EARS 


虽 射 损失 ， 但 同时 会 感应 出 更 大 的 


m? ~ 10?! m3. 


， 在 万 高 达 10 MA 的 聚变 


这 是 


最 大 值 后 快速 衰减 为 零 。 考 虑 到 散 裂 红 
m? 的 缓解 物质 ,未 来 可 利用 SPI 技术 注入 大 量 的 气 


密度 为 


究 了 密度 为 10? m 


FAL (SPD 技术 更 有 可 外 
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是 影响 破裂 后 平台 逃逸 电流 大 小 的 关键 
FELA E. HEARST 
体 通常 有 一 最 优 含量 , FETT A eC EP A A i HTC PAESE TT RS 


或 气氛 混合 气 
|. ERE 
量 降 低 , 不 能 有 效 地 
EE, 虽然 


欧姆 电场 ， 产 生 更 高 的 逃逸 电流 。 
E 入 条 件 下 ， 研 究 发 现 破裂 平台 逃逸 电流 态 受 Tp 


入 气 握 / 氛 混合 气体 的 总 量 也 


\ 握 / 氛 混合 气体 中 的 含量 应 在 0.50% 


量 级 。 在 所 讨论 


~0.70%， 氛 的 注入 量 


T 


级 托 卡 马克 装置 上 ， 注 入 混 


前 MGI 技术 所 不 能 达到 
3 条 件 下 碰撞 
电流 ， 而 是 形成 峰值 逃逸 电流 ， 即 逃逸 电流 随时 间 达 到 


的 。 
耗 散 效 应 对 逃逸 电流 的 影响 。 


实现 破裂 后 向 等 离子 体 注入 


上 高 等 离子 体 电流 


SoM HAR FU 


下 的 聚变 堆 托 卡 马 克 等 离子 体 的 破裂 逃逸 电流 。 


作者 贡献 说 明 HAH: 


beth 


参考 文献 


[1] Smith H, Helander P, Eriksson L G, et al. Runaway electrons and the evolution of the plasma current in tokamak 


究 思 路 、 对 文章 知识 性 


进行 研究 调查 、 数 值 模拟 、 数 据 处 理 、 论 文 初稿 撰写 ; 郑平 卫 : 
内 容 作 批评 性 审阅 、 论 文 指导 与 修改 。 


disruption[J]. Physics of Plasmas, 2006, 13(10): 102502. DOI: 10.1063/1.2358110. [LinkOut ] 


[2] Yoshino R, Nakamura Y, Neyatani Y. Avoidance of VDEs during plasma current quench in JT-60U[J]. Nuclear 
Fusion, 1996, 36(3): 295 — 307. DOI: 10.1088/0029-5515/36/3/i03. [LinkOut ] 

[3] Loarte A, Riccardo V, Martin- Solis J R, et al. Magnetic energy flows during the current quench and termination 
of disruption with runaway current plateau formation in JET and implications for ITER[J]. Nuclear Fusion, 2011, 
51(7): 073004. DOT: 10.1088/0029-5515/51/7/073004. [LinkOut ] 

[4] Li Y L, Sun J Z, Zhang Y P, et al. Simulation of runaway electron generation and diffusion during major 
disruption in the HL-2A tokamak[J]. Fusion Engineering and Design, 2014, 89(7 — 8): 1019 — 1023. DOI: 
10.1016/j.fusengdes.2014.02.074. [LinkOut] 

DOI: 10.1103/physrev.117.329. | LinkOut ] 

[6] Smith H M, Fehér T, Fülóp T, et al. Runaway electron generation in tokamak disruption[J]. Plasma Physics and 
Controlled Fusion, 2009, 51(12): 124008. DOI: 10.1088/0741-3335/51/12/124008. [LinkOut ] 

[7] Fehér T, Smith H M, Fülóp T, et al. Simulation of runaway electron generation during plasma shutdown by 
impurity injection in ITER[J]. Plasma Physics and Controlled Fusion, 2011, 53(3): 035014. DOI: 10.1088/0741 - 
3335/53/3/035014. [LinkOut ] 

[8] Berger E, Pusztai I, Newton S L, et al. Runaway dynamics in reactor-scale spherical tokamak disruption[J]. 
Journal of Plasma Physics, 2022, 88(6): 905880611. DOL: 10.1017/50022377822001209. [LinkOut ] 

[9] Zhang Y P, Tong R H, Yang Z Y, et al. Recent progress on the control and mitigation of runaway electrons and 
disruption prediction in the HL-2A and J-TEXT tokamaks[J]. Reviews of Modern Plasma Physics, 2023, 7(1): 1 — 
39. DOI: 10.1007/541614-022-00110-3. [LinkOut] 

[10] Hesslow L, Embréus O, Vallhagen O, ef al. Influence of massive material injection on avalanche runaway 
generation during tokamak disruption[J]. Nuclear Fusion, 2019, 59(8): 084004. DOI: 10.1088/1741- 
4326/ab26c2. [LinkOut] 

for disruption mitigation experiments[J]. Nuclear Fusion, 2019, 59(6): 066008. DOI: 10.1088/1741-4326/ab136c. 


[LinkOut] 


[12] 计 华 健 , 张 洪 明 , 杨 秀 达 , 等 . EAST 辐射 偏 滤器 运行 模式 下 等 离子 体 杂质 浓度 评估 方法 思 . 核 技术 , 2022, 


45(3):8. DOI: 10.11889/j.0253-3219.2022.hjs.45.030501 


[13] Vallhagen O, Embreus O, Pusztai I, et al. Runaway dynamics in the DT phase of ITER operations in the 
presence of massive material injection[J]. Journal of Plasma Physics, 2020, 86(4): 475860401. DOI: 
10.1017/s0022377820000859, [LinkOut ] 


[14] Hoppe M, Embreus O, Fülóp T. DREAM: a fluid-ki 


ption runaway 
electron simulations[J]. Computer Physics Communications, 2021, 268: 108098. DOI: 10.1016/j.cpc.2021.1 

08098, [LinkOut] 

[15] Breizman B N, Aleynikov P, Hollmann E M, et al. Physics of runaway electrons in tokamaks[J]. Nuclear Fusion, 


2019, 59(8): 083001. DOI: 10.1088/1741-4326/ab1822. | LinkOut ] 


[16] Helander P , Smith H ,T. Fülp,et al. Electron kinetics in a cooling plasma[J].Physics of Plasmas, 2004, 
11(12):5704-5709.DOI:10.1063/1.1812759. 


[17] Papp G, Fülóp T, Fehér T, et al. The effect of ITER-like wall on runaway electron generation in JET[J]. Nuclear 


Fusion, 2013, 53(12): 123017 10.1088/0029-5515/53/12/123017. | LinkOut ] 

[18] Connor J W, Hastie R J. Relativistic limitations on runaway electrons[J]. Nuclear Fusion, 1975, 15(3): 415- 
424. DOI: 10.1088/0029-5515/15/3/007. [LinkOut ] 

[19] Wei Y N, Yan W, Chen Z Y, et al. Dissipation of runaway current by massive gas injection on J-TEXT[J]. 


Plasma Physics and Controlled Fusion, 2020 , 62(2): 025002. DOI: 10.1088/1361-6587/ab52c8. [LinkOut] 


[20] Rosenbluth M N, Putvinski S V. Theory for avalanche of runaway electrons in tokamaks[J]. Nuclear Fusion, 
1997, 37(10): 1355 — 1362. DOI: 10.1088/0029-5515/37/10/03. [LinkOut ] 

[21] Fehér T, Smith H M, Fülóp T , et al. Simulation of runaway electron generation during plasma shutdown by 
impurity injection in ITER[J]. Plasma Physics and Controlled Fusion, 2011, 53(3): 035014. DOI: 10.1088/0741- 
3335/53/3/035014, [LinkOut ] 

[22] Berger E, Pusztai I, Newton S L, et al. Runaway dynamics in reactor-scale spherical tokamak disruptions[J]. 
Journal of Plasma Physics, 2022, 88(6): 905880611. DOI: 10.101 7150022 377822001209. [LinkOut ] 

[23] Summers H P. 2004 The ADAS user manual, version 2.6.[DB/OL] www.adas.ac.uk/manual.php 

[24] Duan X R, Xu M, Zhong W L, et al. Progress of HL-2A experiments and HL-2M program[J]. Nuclear Fusion, 
2022, 62(4): 042020. DOI: 10.1088/1741-4326/ac3be6. [LinkOut] 

[25] Vallhagen O, Pusztai I, Hoppe M, et al. Effect of two-stage shattered pellet injection on tokamak disruption[J]. 


Nuclear Fusion, 2022, 62(11): 112004. DOI: 10.1088/1741 -4326/ac6676. [LinkOut] 


